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Abstract — We propose distributed algorithms to automati- 
cally deploy a group of robotic agents and provide coverage of 
a discretized environment represented by a graph. The classic 
Lloyd approach to coverage optimization involves separate 
centering and partitioning steps and converges to the set of 
centroidal Voronoi partitions. In this work we present a novel 
graph coverage algorithm which achieves better performance 
without this separation while requiring only pairwise "gossip" 
communication between agents. Our new algorithm provably 
converges to an element of the set of pairwise-optimal partitions, 
a subset of the set of centroidal Voronoi partitions. We illustrate 
that this new equilibrium set represents a significant perfor- 
mance improvement through numerical comparisons to existing 
Lloyd-type methods. Finally, we discuss ways to efficiently do 
the necessary computations. 

I. Introduction 

Coordinated networks of robots are already in use for 
environmental monitoring [1] and warehouse logistics [2]. 
In the near future, improvements to the capabilities of 
autonomous robots will enable robotic teams to revolutionize 
transportation and delivery of products to customers, search 
and rescue, and many other applications. All of these tasks 
share a common feature: the robots are asked to provide 
service over a space. The distributed territory partitioning 
problem for robotic networks consists of designing individual 
control and communication laws such that the team will 
divide a space into territories. Typically, partitioning is done 
so to optimize a cost function which measures the quality of 
service provided by the team. Coverage control additionally 
optimizes the positioning of robots inside a territory. 

This paper describes a distributed coverage control algo- 
rithm for a network of robots to optimize the response time of 
the team to service requests in an environment represented 
by a graph. Optimality is defined with reference to a cost 
function which depends on the locations of the agents and 
geodesic distances in the graph. As with all multiagent 
coordination applications, the challenge comes from reducing 
the communication requirements: the proposed algorithm 
requires only "gossip" communication, that is, asynchronous 
and unreliable pairwise communication. 

A broad discussion of partitioning and coverage control is 
presented in [3] which builds on the classic work of Lloyd [4] 
on algorithms for optimal quantizer design through "cen- 
tering and partitioning." The relationship between discrete 
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and continuous coverage control laws based on Euclidean 
distances is studied in [5]. Coverage control and partitioning 
of discrete sets are also related to the literature on the 
facility location or fc-center problem [6]. Coverage control 
algorithms for non-convex environments are discussed in [7], 
[8], [9] while equitable partitioning is studied in [10]. Other 
works considering decentralized methods for coverage con- 
trol include [11] and [12]. 

In [13] the authors have showed how a group of robotic 
agents can optimize the partition of convex bounded subset 
of R d using a Lloyd-type algorithm with pairwise "gossip" 
communication: only one pair of regions is updated at 
each step of the algorithm. This gossip approach to Lloyd 
optimization was extended in [14] to discretized non-convex 
environments more suitable for physical robots. 

There are three main contributions of this paper. First, we 
present a novel gossip coverage algorithm and prove it con- 
verges to an element in the set of pairwise-optimal partitions 
in finite time. This solution set is shown to be a strict subset 
of the set of centroidal Voronoi partitions, meaning the new 
algorithm has fewer local minima than Lloyd-type methods. 
Second, through realistic Player/Stage simulations we show 
that the set of pairwise-optimal partitions avoids many of the 
local minima which can trap Lloyd-type algorithms far from 
the global optimum. Third, we discuss how to efficiently 
compute our new pairwise coverage optimization. 

This paper is organized as follows. Section [TT] defines 
the domain and goal of our algorithm, while [III] contains a 
review of Lloyd-type gossip algorithms from [14]. SectionlTVl 
presents our new algorithm and its properties, while in [VI] 
we detail its computational requirements. In Section [VTll we 
provide numerical results and comparisons to prior methods, 
and we end with concluding remarks in I Villi 

In our notation, R>o denotes the set of non-negative real 
numbers and Z>o the set of non-negative integers. Given a 
set X, \X\ denotes the number of elements in X. Given sets 
X,Y, their difference is X \ Y = {x e X | x £ Y}. A 
set-valued map, denoted by T : X =4 Y, associates to an 
element of X a subset of Y. 

II. Problem formulation 

The problem domain and goal we are considering is the 
same as that in [14]: given a group of N robotic agents with 
limited sensing and communication capabilities, and a dis- 
cretized environment, we want to apportion the environment 
into smaller regions and assign one region to each agent. The 
goal is to optimize the quality of coverage, as measured by 
a cost functional which depends on the current partition and 
the positions of the agents. 



Let the finite set Q be the discretized environment. We 
assume that the elements of Q, which can be thought of 
as locations, are connected by weighted edges. Let G = 
(Q,E,w) be an (undirected) weighted graph with edge set 
E C Q x Q and weight map w : E — >■ K>o; we let w e > 
be the weight of edge e. We assume that G is connected and 
think of the edge weights as distances between locations. 

In any weighted graph G there is a standard notion of 
distance between vertices defined as follows. A path in G 
is an ordered sequence of vertices such that any consecutive 
pair is an edge of G. The weight of a path is the sum of the 
weights of the edges in the path. Given vertices h and k in 
G, the distance between h and k, denoted do(h,k), is the 
weight of the lowest weight path between them, or +00 if 
there is no path. If G is connected, then the distance between 
any two vertices is finite. By convention, da{h,k) = if 
h = k. Note that dc(h, k) = dc(k, h), for any h,k € Q. 

We will be partitioning Q into N connected subsets 
or territories to be covered by individual agents. To do 
so we need to define distances on induced subgraphs of 
G = (Q,E,w). Given I G Q, the subgraph induced 
by the restriction of G to I, denoted G fl I, is the graph 
with vertex set equal to / and with edge set containing all 
weighted edges of G where both vertices belong to /, i.e. 
(Q,E,w)r\I = (QnI,En(I x I),w\ lxl ). The induced 
subgraph is a weighted graph with a notion of distance 
between vertices: given h,k G /, we write di(h,k) := 
dani{h, k). Note that di(h, k) > dc{h, k). 

We then define a connected subset of Q as a subset S C Q 
such that 5^0 and G n S is connected, with C(Q) denoting 
the set of such subsets. We can then define partitions of Q 
into connected sets as follows. 

Definition II.l (Connected partitions) Given the graph 
G = (Q,E,w), we define a connected N— partition of Q 
as a collection p = {pi}fLi of N subsets of Q such that 

« [J?=iPi = Q; 

(ii) p t n pj ; = if i ^ j; 
(in) Pi^lD for all i £{!,..., N}; 
(iv) pi G C(Q) for allie{l,..., N}. 
Let V to be the set of such partitions. 

Remark II.2 When \Q\ > N and G is connected, it is 
always possible to find a connected N— partition of Q. Pick 
N unique vertices from Q to seed the subsets, then iteratively 
grow each subset by adding unclaimed neighboring vertices. 

For our gossip algorithms we need to introduce the notion 
of adjacent subsets. Two distinct connected subsets Pi, pj are 
said to be adjacent if there are two vertices qi, qj belonging, 
respectively, to pi and pj such that (qi,qj) € E. Observe 
that if pi and pj are adjacent then pi Upj G C(Q). Similarly, 
we say that robots i and j are adjacent or are neighbors if 
their subsets pi and pj are adjacent. Accordingly we next 
provide the definition of an adjacency graph, which we will 
use in the analysis of our gossip algorithms. 



Definition II.3 (Adjacency graph) For p G V, we define 
the adjacency graph between regions of partition p as Q{p) = 
({1, . . . , N}, £ (j>)), where G £(p) if Pi and pj are 

adjacent. 

On Q, we define a weight function to be a bounded positive 
function (j> : Q — > M>o which assigns a positive weight to 
each element of Q. Given pi G C(Q), we define the one- 
center function Hi : pi — > R>o as 

H 1 (h; Pi )= ^dp^MMfc)- 

A technical assumption is then needed to define the 
generalized centroid of a connected subset. In what follows, 
we assume that a total order relation, <, is defined on Q: 
hence, we can also denote Q = {1, . . . , \Q\}. 

Definition II.4 (Centroid) Let Q be a totally ordered set, 
and pi G C(Q). We define generalized centroid of pi as 

Cd(pi) := min{argmin'Hi(/i;pj)}. 

hEpi 

In subsequent use we will drop the word "generalized" for 
brevity. Note that with this definition the centroid is well- 
defined, and also that the centroid of a region always belongs 
to the region. With a slight notational abuse, we define Cd : 
V — > Q N as the map which associates to a partition the 
vector of the centroids of its elements. 

With these notions we can introduce the multicenter func- 
tion "Hmuiticentei- : V x Q N -> R> defined by 

N 

^multicenter (p, c) = ^ Hi (Cj] Pj) . 

i=l 

Our motivation for using this cost function is to optimize the 
response time of the robots to a task appearing randomly in Q 
according relative weights <j>. We aim to minimize 'Hmuiticenter 
with respect to both the partitions p and the points c so as 
to minimize the expected distance from this random vertex 
to the centroid of the partition the vertex is in. 

Among the ways of partitioning Q, there is one which is 
worth special attention. Given points c G Q N such that if 
i 7^ j, then Cj 7^ cj, the partition p G V is said to be a 
Voronoi partition of Q generated by c if, for each pi and all 
k G pt, we have c, G Pi and dc(k, Cj) < dc(k, Cj), Vj / i. 

Proposition II.5 (Properties of muiticenter function) Let 

p G V, c G Q , and let p* be a Voronoi partition generated 
by c. Then 

^multicenter (P: Cd(p)) < Hmulticenter 

^multicenter^ i ^ ^multicenter {jPi c). 

These statements motivate the following definition: a par- 
tition p G V is a centroidal Voronoi partition if p is a Voronoi 



partition generated by Cd(p). Based on the multicenter 
function, we define ^expected : V — > R>o by 

^expected (p) = %multicenter(P; Cd(p) ) 

i— 1 a;Gpi 

Observe that ^expected has the following property as an 
immediate consequence of Proposition III. 51 given p G V, if 
p* is a Voronoi partition generated by Cd(p) then 

^expected (p ) 5^ ^expected 

We are now in a position to state the goal of our algorithm: 
solving the optimization problem 

minHexpected(p), 

using only pairwise territory exchanges between agents along 
the edges of £(p). In the literature this optimization is 
typically studied with either centralized control or with 
synchronous and reliable communication between several 
agents, see [3] and [9]. We believe these communication 
requirements are unrealistic for deployed robotic networks, 
and thus are interested in solutions requiring only pairwise 
gossip communication as first explored in [13] and [14]. 

III. Lloyd-type gossip graph coverage 

In this Section we briefly review the discretized Lloyd-type 
gossip coverage algorithm proposed in [14]. 

Lloyd-type Gossip Graph Coverage Algorithm 



At each time t € Z>o, each agent i G {1, . . . , A} maintains 
in memory a connected subset Pi{t). The collection p(0) = 
{pi(0), . . . ,pjv(0)} is an arbitrary connected A— partition 
of Q. At each t G Z>o a communicating pair, say (i, j) G 
£(p(t)), is selected by a deterministic or stochastic process 
to be determined. Assume that i < j. Every agent k g" {i, j} 
sets pk(t + 1) = Pk(t), while i and j perform the following: 
1: agent i transmits its subset Pi(t) to j and vice-versa 
2: both agents compute centroids Cd(pi) and Cd^j), set 
u = pi U pj, and the sets 

Wi^j = {x Epi : d u (x, Cd(pj)) < d u (x, Cdfe))} 
Wj^i = {x G pj : d u (x,Cd(pi)) < d u (x,Cd(pj))} 
W i9 *j = {x G p, U pj : d u (x, Cd(pi)) = d u (x, Cd(pj))} 

3: if Wi^j U Wj^yi = then 

4: Pi (t + 1) := pi(t) and Pj (t + 1) := Pj (t) 

5: else 

6: Pi (t + 1) := ((pi \ Wi^j) U Wj-n) U Wjsu, 

Pj (t + 1) := {{ Pj \ Wj^i) U Wi^j) \ {W m n Pj ) 
7: end if 



Observe that Wi->j (resp. Wj^i) contains the cells of pi 
(resp. pj) which are closer to Gd{pj) (resp. Cd(pj)), whereas 



Wi^j represents the set of the tied cells. In other words, 
when two robots exchange territory using this Lloyd-type 
algorithm, their updated regions {p,(i + l),pj(t + 1)} are 
the Voronoi partition of the set Pi(t) Upj(t) generated by the 
centroids Cd(pj) and Cd(pj), with all tied cells assigned to 
the agent with the lower index. 

Now, for any pair G {1, . . . , A} 2 , i ^ j, we define 
the map : V =$ V as 

T ,, = [P, if (i,3) t S(p) or W^j U W S -» = 
13 \{pi,...,pi,...,Pj,...,p N }, otherwise, 

where % = {{jpi \ W^j) U Wj^i) U Wj^i, and pj = 
((pAW^UW^MW^np,) . 

The dynamical system on the space of partitions is there- 
fore described by, for t G Z>o, 

p(t + l)=T ij (p(t)), for some eS(p(t)), (1) 

together with a rule for which edge is selected at each 
time. We also define the set-valued map T : V =t V by 

T(p) = {T ij ( J3 )\(i,j)e£(j>)}. (2) 

In Theorem ( IIII.21 we summarize the convergence proper- 
ties of the Lloyd-type discretized gossip coverage algorithm. 
To do so, we need the following definition. 

Definition III.l (Uniform and random persistency) Let 

X be a finite set. 

(i) A map a : Z>o — > X is uniformly persistent if 
there exists a duration A G N such that, for each 
x G X, there exists an increasing sequence of times 
{tk}kez> C Z> satisfying t k+1 - t k < A and 
&{tk) = x for all k G 1>q- 

(ii) A stochastic process a : Z>o — > X is randomly 
persistent if there exists a probability p G ]0, 1[ such 
that, for each x G X and for each t G Z>o 

F[a(t+l) = x\a(t),...,a(l)] > p. 

Theorem III.2 (Convergence under persistent gossip) 

Consider the Lloyd-type discretized gossip coverage 
algorithm T and the evolutions p : Z>o — > V defined by 

p(t+l)=T a(t) (p(t)), forteZ> Q , 

where a : Z>o — > {(*,.?) G {1, . . . , A} 2 | i ^ j} is either a 
deterministic map or a stochastic process. Then the following 
statements hold: 

(i) ;/ a is a uniformly persistent map, then any evolution 
p converges in a finite number of steps to a centroidal 
Voronoi partition; and 

(ii) if a is a randomly persistent stochastic process, then 
any evolution p converges almost surely in a finite 
number of steps to a centroidal Voronoi partition. 

Convergence to a centroidal Voronoi partition puts this 
pairwise optimization algorithm in the same class with 
centralized Lloyd methods. However, there can be a great 
number of possible centroidal Voronoi partitions for a given 



discretized environment and, in our experience, the algorithm 
too frequently converges to a suboptimal solution. This issue 
motivates the developments in the following Section. 

IV. PAIRWISE- OPTIMAL GOSSIP COVERAGE ALGORITHM 

In this section we present a novel gossip graph coverage 
algorithm which improves on the performance achievable by 
the Lloyd-type algorithm reviewed in the previous Section. 
First we introduce the following notational definition. 

Given p £ V and distinct components pi and pj, let P Pi up 
denote the set of all distinct pairs of vertices in pi Lipj. We 
assume that the elements in P Pi u Pj are sorted lexicographi- 



cally. In formal terms, if PiUpj = {hi,. 



.h 



biU Pj | 



} where 



h s < h s+ i for s = {1, . . . , \p.i Upj\- 1}, then 
Pptupj = {(hi,h 2 ), (hi, h 3 ), . . . , (hi, h\ P{Upj \), (h 2 ,h 3 ), . . 



■ j (h\ Pi Upj\-l: fypiUpj|)} • 



Notice that \P„. 



_ \ Pi U Pj \ (\PiUPj\-l) 
PiUpj I — 2 

The pairwise-optimal discretized gossip coverage algo- 
rithm is defined as follows. 



Gossip Optimal Graph Coverage Algorithm 

At each time t £ Z>o, each agent i £ {1, . . . , N} maintains 
in memory a connected subset pi(t). The collection p(0) = 
{pi(Q), . . . ,Pn(0)} is an arbitrary connected N— partition 
of Q. At each t £ Z>o a communicating pair, say (i,j) £ 
£(p(t)), is selected by a deterministic or stochastic process 
to be determined. Assume that i < j. Every agent k £ {i,j} 
sets pk(t + 1) = Pk(t), while i and j perform the following: 
1: agent i transmits its subset Pi(t) to j and vice-versa 
2: both agents compute the set u = pi Upj, the function 
V : P u -> M >0 



V((a,b)) 



{d u (h,a),d u (h,b)} . 



and find the set St> = argmirw ^gp^ {T>(h s , h r )} 
3: assuming St> is lexicographically ordered, both agents 

pick the first pair in S-p, say (a*, b*) 
4: both agents compute the sets 

W a * = {x £ u : d u (x, a*) < d u (x, b*)} 
Wb* = {x £ u : d u (x,b*) < d u (x,a*)} 



if Ui(Cd(W a .);W a . 
miCdip^t));^)) 
ft(t+l) = W a ; 
else 

Pi (t + l)= Pi (t), 
end if 



+ fti(Cd(w 6 .);W6*) < 

-Hi(Cd( Pj (t));j3 3 (t)) then 
Pj (t + l) = W b . 

Pj (t + l)= Pj (t) 



Remark IV. 1 The lexicographic rule for picking (a*,b*) 
used here makes the dynamical system represented by the 
algorithm deterministic and well-defined. In practice any pair 
in St> can be chosen. 



The following Proposition characterizes a desirable prop- 
erty of the regions W a * and Wb* ■ 

Proposition IV.2 The regions W a * and Wb- as defined in 
Gossip Optimal Graph Coverage Algorithm, are connected. 



Thanks to Proposition IIV.2I we have that, for any pair 
(i,j) £ {1, • • . ,N} 2 , i 7^ j, we can define the map : 
V^Vby 



T ij(p) 



if (i,j)££(p) 
Si, ■ ■ ■ ,Pj,- ■ ■ ,Pn}, 



otherwise, 



where pj = W a * and pj = Wb* ■ Therefore, the dynamical 
system on the space of partitions can again be described in 
compact form by (Q]i and (f2]i but with this new T^ . 

To establish the convergence properties of the pairwise- 
optimal discretized gossip coverage algorithm we need the 
following definition. 

Definition IV.3 A partition p £V is said to be a pairwise- 
optimal partition if, for every (i,j) £ £(p) 

Wi(Cd( Pl );p i )+«i(Cd(p i );Pi) = 

min < S~] min{d PiUp (a, k),d PtUp Ab, k)} 

a.b£piUVj * — ' 

{ke Pl u Pj 

The following Proposition characterizes the set of the 
pairwise-optimal partitions. 

Proposition IV.4 Let p £ V be a pairwise-optimal partition. 
Then p is also a centroidal Voronoi partition. 

We are now ready to state the main result of this paper. 

Theorem IV.5 (Convergence under persistent gossip) 

Consider the pairwise-optimal discretized gossip coverage 
algorithm T and the evolutions p : Z>o — > V defined by 

p(t + l)=T a{t) (p(t)), fort£Z> , 

where a : Z>o — > {(i,j) £ {1, ■ • ■ , N} 2 | i ^ j} is either a 
deterministic map or a stochastic process. Then the following 
statements hold: 

(i) if a is a uniformly persistent map, then any evolution 
p converges in a finite number of steps to a pairwise- 
optimal partition; and 

(ii) if a is a randomly persistent stochastic process, then 
any evolution p converges almost surely in a finite 
number of steps to a pairwise-optimal partition. 

In contrast to the Lloyd-type algorithm in Section [III] this 
new pairwise-optimal algorithm converges to an element in 
the set of pairwise-optimal partitions. This result is a useful 
improvement as the set of pairwise-optimal partitions is a 
subset of the set of centroidal Voronoi partitions, so that 
the new algorithm can avoid many of the local minima in 
which Lloyd-type methods get stuck. We will demonstrate 
this improvement in Section IVIII 



V. Convergence proofs 



we have that 



In this Section we provide proofs of the results stated in 
the previous Section, starting with Proposition IIV.2I 

Proof: [Proposition II V.2I To prove the statement of the 
Proposition we show that, given x G Wb* , any shortest path 
in pi Upj connecting x to b* completely belongs to Wb* ■ We 
proceed by contradiction. Let s x ,b* denote a shortest path in 
Pi Upj connecting x to b* and let us assume that there exists 
m G s x .b* such that m G W a * ■ For m to be in b* means that 
d PiUpj (to, b*) < d PiUpj {m,a*). This implies that 

d Pl u Pj (x, a*) = d PtUpj (to, a*) + d PtUpj (x, to) 

> dpiUpj ( m , b *) + d Pl u P] (x, m) 

> d PiUpj (x,b*). 

This is a contradiction for x G b*. Similar considerations 
hold for W a *. ■ 
To prove Proposition IIV.4I we need to review the notion 
of a centroidal Voronoi in pairs partition [14]. 

Definition V.l Consider the connected graph G = 
(Q,E,w).. Let p G V and let Q{p) be the associated 
adjacency graph. Then p G V is said to be centroidal Voronoi 
in pairs if, for any i G {1, . . . , N} we have that 

d Pt u P] (k, Cd(pi)) < d PtUP] (k, Cd(pj)) (3) 

for all k € pi and for all j such that G Sip). 

The following Lemma [14] states a useful equivalence 
property between the set of centroidal Voronoi in pairs 
partitions and the set of centroidal Voronoi partitions. 

Lemma V.2 A partition p G V is centroidal Voronoi in pairs 
if and only if it is centroidal Voronoi. 



We are now ready to prove Proposition IIV.4I 

Proof: [Proposition lIV.4l We will show that if a partition 
p G V is a pairwise-optimal partition then p is also centroidal 
Voronoi in pairs. Then from Lemma IV.2I it will follow that 
p is also centroidal Voronoi. 

We proceed by contradiction. Let p G V be a pairwise- 
optimal partition, let pi and pj be two adjacent regions, and 
let us assume that there exists x G pi such that 

d Pi \j P j i%> Cd(pi)) > d PtUpj (x, Cd(pj)). (4) 

Let s Xt cd( P ) denote a shortest path in pi U pj connecting 
x to Cd(pj). Without loss of generality let us assume that 
s x ,cd( P] ) = (x,yi,V2, ...,y m ,Cd(pj)), for a suitable m- 
tuple in pi U pj. Now let s^ca^) denote the subset of 
s x,cd( Pj ) consisting of only the vertices belonging to pi, i.e., 

s i l ,cdfo) = {y e s x,cd( Pj )\y epi}. 

Observe that from d3]i it follows that 



d Pt u Pj (y, Cd{p l )) > d PiUP] (y, Cd(pj)) 



for all y G s K 



(0 



Gd( Pj ) 



Therefore, by letting 



HiiCdip^Pi) + u 1 (cd{p j y ) p j ) 

K-H^Cdip^Pi) + Hi(Cd( Pj );p j ) 
<Ki(Cd(p i );p i )+^ 1 (Cd(p i ); % ). 

This last inequality contradicts the fact that p is a pairwise- 
optimal partition. ■ 
To provide the proof of Theorem IIV.5I we need the 
following intermediate result. 

Lemma V.3 Let T : V =4 V be the point-to-set map defined 
for the Gossip Optimal Graph Coverage Algorithm. If pit) G 
V and p(t + 1) G T{p(t)) \ p(t), then H ewected (p(t + 1)) < 

Proof: Without loss of generality assume that 
is the pair selected at time t. For simplicity let us denote 
p(t + 1) and pit), respectively, by p + and p. Then 



^expected iv ) ^expected ip) 



);pt)+Ki(Cd(p+);pt) 



- iHiiCdip^pi) + H i (Cd(p j y iPj )). 

According to the definition of the Gossip Optimal Graph 
Coverage Algorithm we have that if pf ^ Pi, P^ ^ Pj, then 

Wi(Cd(p+);pt)+%(Cd( P t);p+) 

<^ 1 (Cd( P i);Pi)+^i(Cd(p j );p i ) 

from which the thesis follows. ■ 
Now we can prove the main result of this paper. 

Proof: [Theorem IIV.5I The proof of point (i) of the 
Theorem is based on on verifying the assumptions (i), (ii), 
(iii), (iv) of Theorem 4.3 from [13]. 

We start with the following topological consideration. Let 
A, B G Q, AAB be their symmetric difference, and let \A\ 
be the "points counting measure" of A, that is the number 
of elements of A. Then if we define d&.iA,B) = \AAB\ , 
we have that g?a is a distance on the set of the subsets of 
Q, which we denote 2^. Since the points counting measure 
takes integer values, 2^ a discrete topological space. Next, 
consider the product space (2®) N and let q = {qi, . . . , q^} 
and q = {<7i, • • • , 5jv} be two its elements. Then, thanks 
to results valid for product topological spaces, we have that 



the function d^'iq^) = YliLi d (?i> Qi) ' s a distance on 
i2®) N . Hence, (2®) N is also a discrete topological space. 
Observe that P is a subset of (2 ( 9) iV . Since 2®, [2^) N 
and V are finite, they are also trivially compact. Moreover, 
since the algorithm T : V =4 V is well-defined, we have 
that V is strongly positively invariant. This fact implies that 
assumption (i) of Theorem 4.3 from [13] is satisfied. 

From Lemma I V.3 1 it follows that assumption (ii) is also 
satisfied, with ^expected playing the role of the function U. 

Next, the fact that all maps from a discrete space are 
continuous implies the continuity of the maps and 
"^expected, thus guaranteeing the validity of assumption (iii). 
Assumption (iv) holds true for the hypothesis made in the 



Pi 



Pi \ S x,Cd( Pj ) 



and 



(i) 



Pj =PiUs;' c% ) 



Fig. 1. At left, ten robots are positioned at the centroids of their initial partitions in a non-convex environment. The boundary of each agent's partition is 
drawn in a different color. The middle shows the final solution of a centralized Lloyd optimization with a total cost of 624.1m. At right is a final solution 
for the pairwise-optimal discretized gossip coverage algorithm starting from the final centralized Lloyd solution, with an improved cost of 610.3m. 



statement of Theorem IIV.5I that a is a uniformly persistent 
map. Hence we are in the position to apply Theorem 4.3 
from [13], and conclude convergence to the intersection 
of the equilibria of the maps Xy, which, according to the 
definition of 2y coincides with the set of the pairwise- 
optimal partitions. Hence, since that set is finite, we can 
argue that the system converges in finite time to one pairwise- 
optimal partition. 

A proof of (ii) in Theorem II V. 5 1 can be made in a similar 
way using Theorem 4.5 from [13] inplace of 4.3. ■ 

VI. Scalability Properties 
In this Section we discuss how to compute the optimal 
new centroids (a*,b*) £ -P Pi up f° r the Gossip Optimal 
Graph Coverage Algorithm, as well as a sampling method 
to greatly reduce the computational complexity. Determining 
(a* , b*) requires an exhaustive search over all possible pairs 
of vertices for a pair with the lowest cost to cover pi U Pj . 
Using Johnson's all pairs shortest paths algorithm as a first 
step to compute a pairwise distance matrix for PiUpj, (a*, b*) 
can be found in 0(\pi\ 3 ) with a memory requirement of 
0(\pi\ 2 ) (as pi Upj is of order \p t \). 

For mobile robots with limited onboard resources, these 
requirements may be too steep. In such circumstances, we 
propose instead to sample pairs of potential new centroids 
when applying the algorithm. The first sample pair would 
be the agents' previous centroids, with the rest drawn at 
random from the set P Pi u Py For m sample pairs, this 
approach requires 2m calls to Dijkstra's one-to-all shortest 
paths algorithm to find the pair with the lowest cost, for 
a total time complexity of 0(m\pi \ \og(\pi\)) and a memory 
requirement of 0([pj|). It is also worth noting that if all edge 
weights in G are equal (i.e., for a regular grid discretization), 
then a breadth-first-search approach can replace Dijkstra and 
the time complexity drops to 0(m\pi\). This sample-based 
approach greatly reduces the time and memory requirements, 
and will still converge to a pairwise-optimal partition. 

VII. Simulation Results 

To demonstrate the utility of the proposed gossip coverage 
algorithm, we implemented it and other coverage algorithms 



using the open-source Player/Stage robot software system 
[15] and the Boost Graph Library (BGL) [16]. All results 
presented here were generated using Player 2.1.1, Stage 
2.1.0, and BGL 1.34.1. A non-convex environment was 
specified with a bitmap and overlaid with a 0.1m resolution 
occupancy grid, producing a lattice-like graph with all edge 
weights equal to 0.1m. To compute distances on graphs with 
uniform edge weights we extended the BGL implementation 
of breadth-first-search with a distance recorder event visitor. 

Figure [T] provides one example of how convergence to the 
set of pairwise-optimal partitions represents an improvement 
over Lloyd-type methods. On the left is the non-convex 
environment used for all results in this Section, as well as 
an initial partition for 10 robots. The middle panel shows 
the result of a centralized Lloyd optimization where all 
agents iteratively update their centroid and partition based 
on the centroids of all other agents. The final equilibrium 
for this centralized method is a centroidal Voronoi partition 
with a total cost of 624.1m. However, this solution is not a 
pairwise-optimal partition. We started the pairwise-optimal 
discretized gossip coverage algorithm using this partition as 
the initial condition. After 70 pairwise exchanges between 
random agent pairs, the algorithm reached the pairwise- 
optimal partition at right which has a lower cost of 610.3. 
The best known 10-partition of this map has cost 610.0. 

Figure [2] compares 100 simulations of the Lloyd-type 
and pairwise-optimal discretized gossip coverage algorithms. 
We started each simulation from the initial partition in 
Fig. 03 and used different random sequences of agent pairs 
for each case. The new pairwise-optimal algorithm shows 
marked improvement over the Lloyd-type gossip algorithm, 
achieving a total cost within 2% of the best known cost in 
85% of trials, while also avoiding the worst local minima. 

Figure [3] compares final cost histograms for 10 random 
initial conditions for the same environment. Each initial 
condition was created by selecting unique starting locations 
for the agents uniformly at random, and using these locations 
to generate a Voronoi partition. The histograms compare 
20 simulations of Lloyd-type and the new pairwise-optimal 
gossip coverage starting from the same initial partition but 
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Fig. 2. Histograms of final total costs for 100 simulations of pairwise- 
optimal (gray) and Lloyd-type (black) gossip coverage. All runs started from 
the initial partition in Fig. [T] with different random sequences of agent pairs. 
The dashed red line shows the final cost using a centralized Lloyd algorithm. 

with different random orders of exchanges. The centralized 
Lloyd final cost for each initial condition is also shown. 
The new pairwise-optimal method outperforms both Lloyd- 
type methods for all 10 tests, although the Lloyd-type gossip 
method is close in two of the trials. These results illustrate 
that convergence to a pairwise-optimal partition represents 
a significant performance enhancement over classic Lloyd 
methods. Interestingly, the Lloyd-type gossip algorithm also 
substantially outperforms the centralized version in 8 of the 
trials. We speculate this is due to the gossip method taking 
trajectories through the space of connected TV— partitions 
which are not possible for the centralized approach. 

VIII. Conclusions 

We have presented a novel distributed coverage control 
algorithm which requires only pairwise communication be- 
tween agents. The classic Lloyd approach to optimizing 
quantizer placement involves iteration of separate centering 
and Voronoi partitioning steps. In the graph coverage domain, 
the separation of centroid and partition optimization seems 
unnecessary computationally for gossip algorithms. More 
importantly, we have shown that improved performance can 
be achieved without this separation. Our new pairwise- 
optimal discretized gossip coverage algorithm provably con- 
verges to a subset of the centroidal Voronoi partitions which 
we labeled the set of pairwise-optimal partitions. Through 
numerical comparisons we demonstrated that this new subset 
of solutions avoids many of the local minima in which both 
gossip and centralized Lloyd-type algorithms can get stuck. 
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